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We discuss simple models which predict the existence of significant gamma-ray fluxes from dark 
matter annihilation. In this context the dark matter candidate is a Majorana fermion with velocity- 
suppressed tree-level annihilation into Standard Model fermions but unsuppressed annihilation into 
photons. These gamma lines can easily be distinguished from the continuum and provide a possibility 
to test these models. 


I. INTRODUCTION 

The nature of the dark matter (DM) component of 
our Universe is a long standing mystery and one of the 
most challenging questions in fundamental physics. It 
has been argued for a long time that the observation of 
monochromatic gamma-ray lines from dark matter anni¬ 
hilation would be a “smoking gun” signature for the par¬ 
ticle nature of dark matter. This statement is correct in 
the sense that it is unlikely that this signature could have 
an astrophysical origin. However, the question whether 
in any DM model such a feature is generic is much more 
subtle. 

In order to clearly identify a gamma line from dark 
matter annihilation the contribution of the dark matter 
annihilation to final state radiation (FSR) and any an¬ 
nihilation channel contributing to the gamma continuum 
have to be small and the annihilation into one or more 
of the final states 77, /17 and Z7 has to be large. Typ¬ 
ically, it is difficult to have such a scenario because the 
final state radiation occurs at tree level while the gamma 
lines are possible only at the quantum level, see Fig. 1. 
For a review on DM annihilation into gamma rays see 
Ref. [1], 

In the simplest dark matter models where the dark 
matter relic density is defined by the annihilation into the 
Standard Model (SM) fermions one generically expects 
to have a large contribution to final state radiation. For 
example, if the dark matter candidate has spin zero and 




Figure 1. DM annihilation with photons in the final state. 


annihilates through the Higgs into SM fermions one can 
explain the relic density through thermal freeze-out but 
at the same time the contribution to final state radiation 
tends to be large. In this case one can find regions in 
the parameter space where the FSR is suppressed 1 but 
generically it is difficult to see a gamma line. 

The case of Dirac dark matter is similar to the scalar 
dark matter scenario. For example, if the Dirac DM can¬ 
didate is charged under an Abelian gauge force one can 
explain the relic density through the interactions with 
the new gauge boson. Typically, the tree level annihila¬ 
tion is not suppressed and therefore one expects to have 
a large contribution to final state radiation. 

In the context of the Minimal Supersymmetric Stan¬ 
dard Model the lightest neutralino, which is a Majorana 
fermion, can be a cold dark matter candidate [5] . A case 
with interesting 7-ray phenomenology is when the neu¬ 
tralino is bino-like and the annihilation through the t- 
channel is velocity suppressed. However, in this set-up 
virtual internal bremsstrahlung (VIB) from the mediat¬ 
ing squarks and sleptons leads to a gamma ray contin¬ 
uum. This contribution spoils the gamma line visibility. 
It has been argued in the literature that the VIB process 
itself can mimic a line signal, but it is only possible in a 
strongly tuned mass degenerate scenario. See Refs. [2, 6- 
14] for recent studies of gamma spectra in different dark 
matter models. 

A case where one could easily have a visible gamma line 
corresponds to the scenario where the tree level annihila¬ 
tion of a Majorana fermion into SM fermions is velocity 
suppressed. Since the dark matter velocity today is very 
small one could hope that the contribution to final state 
radiation is small but the rate for gamma lines is large. 
In this way the velocity suppression compensates for the 
loop factor present in the gamma lines. 

In this letter we investigate a simple scenario where 
the dark matter candidate is a Majorana fermion which 
naturally gives rise to a significant gamma line and neg¬ 
ligible contribution to the continuum. In this way one 
could hope that experiments such as H.E.S.S. looking for 
gamma lines from dark matter annihilation could cut into 
the parameter space of this model. 


1 For examples, see Refs. [2—4]. 
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II. MAJORANA DARK MATTER 

Let us consider a model where the dark matter is a 
Majorana fermion % charged under a new local gauge 
symmetry. Once the new gauge boson is integrated out 
one could have the effective interaction 

£dm D ^X7 / Vx/lV7 5 /- (1) 

Here / could be a SM fermion or a new fermion needed for 
anomaly cancellation in models with extra gauge symme¬ 
tries. Once we dress this operator to study the annihila¬ 
tion XX ~± 77; one finds that this channel is not velocity 
suppressed. The annihilation \X ^7 is suppressed 
when / is a new fermion, because the Yukawa couplings 
for the new fermions inside the loop have to be small. In 
Ref. [15] it was shown that new chiral fermions with large 
Yukawa couplings to the SM Higgs significantly change 
the h —> 77 cross section and are therefore ruled out. 
Thus, new fermions needed for anomaly cancellation need 
to be vector-like under the SM gauge group and obtain 
the dominant part of their mass from a new Higgs field, 
not from the SM Higgs. The case of \X —> Zj is similar 
to the annihilation into gamma gamma. 

The Majorana fermion could have only the following 
interaction 

£dm D ^X7 / "7 5 X/7m/- (2) 

In this case the annihilation into 77 is not present and 
h'y is velocity suppressed, while the annihilation into Z 7 
is not velocity suppressed. See Table I for the summary 
of the features of all possible channels leading to gamma 
lines. In general one could have a model where both 
interactions in Eqs. (1) and (2) are present and one could 
hope to observe the 77 and Z 7 channels. In what follows, 
we discuss a model where the dark matter annihilation 
into gamma lines is naturally large and could provide a 
possibility to test this model. 

III. BARYONIC DARK MATTER 

We investigate the visibility of the gamma lines from 
dark matter annihilation in models where the baryon 
number is a local gauge symmetry spontaneously bro¬ 
ken at the low scale. These models were proposed and 
studied in Refs. [16-20] and are based on the gauge group 


In the two realistic models [16, 17] one has new fermions 
with baryon number which are necessary to define an 
anomaly-free theory. In this context the lightest field 
with baryon number in the new sector is neutral and can 
describe the cold dark matter in the Universe. In both 
cases the dark matter candidate \ can be a Majorana 
fermion if the baryon number for the new fields is ±3/2. 

The interaction of the dark matter candidate with the 
leptophobic gauge boson Zb is given by 

-CD ^5bX7m7 5 X^b; (4) 

where g B is the new gauge coupling. Here we neglect for 
simplicity the kinetic mixing between SM hypercharge 
and baryon number. Notice that all the SM quarks will 
feel this new force associated with the baryon number, 

-CD ^g n qYqZ^. (5) 

The new fermions /' needed for anomaly cancellation ac¬ 
quire mass from a new Higgs and will have axial couplings 
with the Zb gauge boson 

-CD±\gB]’l^f'Z» B . ( 6 ) 

Notice that once the Z B is integrated out one finds the in¬ 
teraction in Eq. (1) with / being the new fermions needed 
for anomaly cancellation. 

In order to investigate the predictions for the gamma 
lines we will focus on the model proposed in Ref. [16]. In 
this context one has the additional fields 

~ (1, 2,1/2, -3/2), ~ (1, 2,1/2,3/2), 

VR ~ (1; 1, -1, -3/2), and g L ~ (1,1,-1, 3/2), 

with the transformation properties under Gb given. 
These new fields together with the dark matter candi¬ 
date acquire mass once the scalar field S B ~ ( 1 , 1 , 0 , —3) 
obtains a vacuum expectation value ( S B ) = v B /V 2 and 
breaks the local baryon symmetry. The v B also deter¬ 
mines the mass of the Z B as 

M Zb = 3g B VB- (7) 

The relevant interactions of the new fermion fields are 
given by 


G b = SU(3) C ® SU{2) l ® U(1) Y 0 U{l) B - (3) 


Table I. Effective operators for Majorana DM and expected 
strength of the possible gamma lines. 

Operator 77 Z 7 /17 

X7 M 7 5 X/7e7 5 / OK OK suppressed 
X7 m 7 5 X/7(j/ OK suppressed 


— C D A^H/lT^Sb + XtiVrUlSb + A xXrXlSb 

+ y XlXlSb + y XrXrS b + h - c - ( g ) 

Once the U(l) B symmetry is broken one has two neutral 
Majorana fields in the xl~XR sector. The mass matrix 
in the (xl, Xl) basis is 



Ai A x \ 
A x A 2 J 


( 9 ) 
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M x [GeV] 

Figure 2. Relic density vs. the dark matter mass for three 
choices of parameters. Note that the points in parameter 
space where the correct relic density is achieved correspond 
to the benchmark scenarios I III. The solid light-blue line rep¬ 
resents the currently allowed range for the DM relic density, 
Q DM /i 2 = 0.1199 ± 0.0027 [22], 


Neglecting the annihilation through the Higgs fields, in 
Fig. 2 we show the relic density vs. the dark matter mass 
for the benchmark scenarios in Table II. The calculation 
of the relic density is performed in the usual way [ 21 ]. 
The correct relic density can be achieved far from the 
resonance even if the annihilation cross section is veloc¬ 
ity suppressed because the DM velocity is about 0.3c at 
freeze out and the suppression is therefore small. 

In this model the predictions for direct detection de¬ 
pend on the mixing angle 9 between the SM Higgs and the 
new Higgs breaking the local baryon number, since the 
direct detection cross section through the Zb is velocity 
suppressed; see Ref. [20] for more details. As it has been 
shown there, one can have large rates when the mixing 
angle is close to the experimental bound, 9 < 0.35. We 
have checked the possible bounds from the Large Hadron 
Collider (LHC) searches for this model and would like to 
mention that the results shown in this paper are in agree¬ 
ment with the current bounds. See for example Ref. [23] 
for a detailed discussion of the LHC bounds. 

IV. VISIBLE GAMMA LINES 


Given the parameter constellation with the following hi¬ 
erarchy A 2 A x Ai we have a seesaw in the dark 

A 2 

matter sector leading to mi « and m 2 ~ ^ 2 ^- 

Thus Xi is the lightest particle in this sector, is stable 
and can be the cold dark matter candidate. Note that 
the stability of the dark matter is due to a remnant 
symmetry in the new sector. Furthermore, in this set-up 
X 2 is significantly heavier than xi and co-annihilations 
do not play any significant role in the freeze-out process. 
Additionally no co-annihilations with the other neutral 
fermionic states occur, as we work in the regime where 
the Yukawa couplings of the new fermions are small and 
thus there is practically no mixing with the other neutral 
states. 

In our study we focus on the scenario where the dark 
matter is the lightest new field in the theory. We show the 
predictions assuming that the relic density can be ther¬ 
mally produced through the freeze-out of our dark matter 
candidate. The relevant annihilation channel then is the 
annihilation to SM quarks, XX ~A %b “A 99i an d the 
corresponding cross section is given by 


c (XX -A qq) = 


4 yjs - 4M*y/s - 4M2 ( s + 2 Ml) 


9b 

167TS 


( s _ M |j + r| B Mf B 


( 10 ) 

where s is the square of the center-of-mass energy. The 
s-wave contribution to this annihilation process is van¬ 
ishing identically, such that it is p-wave dominated and 
hence velocity suppressed. The other possible annihila¬ 
tion channels through the two Higgs bosons have a double 
suppression: they are suppressed by the mixing angle be¬ 
tween the two Higgs fields and also velocity suppressed. 


Recently, we have discussed the dark matter annihi¬ 
lation into gamma lines in simplified dark matter mod¬ 
els [13]. The calculations performed in Ref. [13] are gen¬ 
eral and can be used for the baryonic dark matter model 
discussed here. In this model, there are three possible 
gamma lines, 

XX -A Z*B -T 77 , hj, Zj, (11) 

and the effective couplings Zb'J'Y, Zsh'y, and ZbZ 7 
are generated by loops of SM quarks and/or the heavy 
fermions T and 17 . For relic density and for DM direct de¬ 
tection, the model has three free parameters: M x , Mz B , 
and gs- For the indirect detection, one has two extra 
parameters in the model, My and M^, the masses of the 
heavy new fermions, which gives a total of five free pa¬ 
rameters. We define three benchmark scenarios for this 
set of parameters in Table II, which we will use through¬ 
out this letter. It is important to remind the reader that 
in the case of the annihilation into two gammas the en¬ 
ergy of the gamma line is equal to the dark matter mass. 
In the case of the annihilation into h'y and Z 7 the energy 
of the gamma line is 

E ->=( [ - §i) ■ < i2 > 


Table II. Benchmark scenarios. 


Scenario M x [GeV] M Zb [GeV] g B M v [GeV] My [GeV] 


I 

300 

360 

0.445 

300 

600 

II 

500 

600 

0.580 

1000 

1500 

III 

420 

800 

0.240 

420 

500 




















4 


Table III. Cross sections in 10 30 cm 3 s 1 . 


Scenario (<tv) 77 

{<jv) Zi 

(vv)h~/ 

(™)bb 

I 

11.5 

1.1 

8.0 x 10 -7 

0.023 

II 

0.021 

15.2 

6.1 x 10 -7 

0.024 

III 

0.047 

0.001 

1.3 x 10~ 6 

0.048 


with X = h or X = Z for the corresponding channel. 

The cross section for the annihilation to two gammas 
is given by 


velocity for the final state radiation will be of the form 

(<™)fsr = A —| + Bv 2 + 0(v A ), (15) 

M x 

with some coefficients A and B. Notice that the first 
term is suppressed by the quark masses, and the second 
term by the small DM velocity today. 

The differential gamma flux from the annihilation into 
two photons is given by 


d<f > 7 

dE 1 


^ann (0^)77 
47 tM£ 


G(E^/w,M x ), 


(16) 


c(xx 77) 


X 


81a 2 9%-^x 

" 327T 3 M| B (m| b +t| b ) 

[MlCo{s-,Mv) - M%Co(s-,M v )] 2 



(13) 


where a = e 2 /4n. The scalar Passarino-Veltman func¬ 
tion Cq is given by 


C 0 (s;m) = — In 2 
2 s 



(14) 


We have used Package-X [24] to calculate the one-loop 
integrals. Some comments regarding the dark matter an¬ 
nihilation into gamma gamma are in order. The cross 
section is not velocity suppressed, and therefore the rate 
can be large. The effective coupling Zb 77 relies on a 
non-vanishing axial coupling of the fermions running in 
the loop to the Zb, and therefore vanishes for the SM 
quarks. This means that only the new fermions needed 
for anomaly cancellation define the annihilation into pho¬ 
tons. 

In Table III, we give the annihilation cross section to 
two gammas for the three benchmark scenarios defined in 
Table II, together with the other annihilation cross sec¬ 
tions for indirect detection. The prediction of the cross 
section for the annihilation into gamma gamma is given 
in Fig. 3, together with the experimental bounds from 
Fermi-LAT [25] and H.E.S.S. [26]. As one can appre¬ 
ciate, one can achieve a rate close to the experimental 
limits when the new fermions are light. The dark matter 
annihilation into a gamma and the SM Higgs, XX ^75 
is velocity suppressed. See Ref. [13] for the details. 


V. GAMMA SPECTRUM 

The key feature of this type of model that allows for a 
visible gamma line is the fact that the DM annihilation 
into two photons is unsuppressed, whereas the tree-level 
annihilation to SM fermions and the final state radia¬ 
tion are suppressed. Generically the cross section times 


where we take the gamma line as a delta peak and model 
the detector resolution by a Gaussian 

1 (E^-Eq) 2 

G(E 7 ,£/u,,E 0 ) = e~^^. (17) 

V2ttE 0 (^/uj) 

Here the parameter £ is a measure of the energy resolu¬ 
tion of the detector and will vary typically between 1 % 
and 10%, and the factor ui = 2-^2 log 2 ss 2.35 determines 
the full width at half maximum, with the standard de¬ 
viation given by <to = Eq^/co. The J-factor J ann is the 
integral over the DM profile, which contains the astro- 
physical assumptions about the DM distribution in the 
galaxy and thus all astrophysical uncertainties. 

In Fig. 4 we show the differential spectrum for the dark 
matter annihilation into gamma rays in the benchmark 
scenario I in Table II. The dominant contribution to the 
continuum in this scenario is FSR from the annihilation 


fi DM h 2 = 0.1199 ± 0.0027, M Zb = 1.2 M x 



Figure 3. Cross section for the DM annihilation into two gam¬ 
mas. We show two setups which have Mz B = 1.2 M x . Setup 1 
(black solid line) is defined by the parameter choice = M x , 
My = 2 M x . Setup 2 (green dashed line) is defined by the pa¬ 
rameter choice M v = 2 M x , M* = 3 M x . The gauge coupling 
gB is chosen such that in each point the dark matter relic 
density has the current value, Ddm/i 2 = 0.1199 ± 0.0027 [22]. 
Experimental limits on the cross section from Fermi-LAT [25] 
and H.E.S.S. [26] are given. 
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Figure 4. Differential spectrum of the annihilation of DM 
into gamma rays for benchmark scenario I. The three curves 
show the spectra with experimental resolutions of 10 % (blue 
dotted), 5% (red dashed) and 1% (green solid) in energy. We 
use the J-factor for the R3 region-of-interest by the Fermi- 
LAT collaboration, J ann = 13.9xl0 22 GeV 2 cm “ 5 [27]. The 77 
and Z 7 lines are included. The continuum spectrum consists 
of the final state radiation, since the other channels leading 
to a continuum gamma spectrum are velocity suppressed. 

into b quarks. Since the maximal energy of the FSR 
photons is given by 

^FSR.max = M x (l - M\/ M 1 2 3 4 5 6 ) , (18) 

the endpoint of the FSR gamma spectrum from the top 
quark will be far from the gamma line unless the DM 
is much heavier than the top quark. Then, however, 
the FSR cross section, see Eq. (15), will be helicity sup¬ 


pressed by M t 2 /M 2 -C 1. 

Notice that the parameter choice in Fig. 4 is repre¬ 
sented by the black line in Fig. 3 where the predicted 
cross section is close to the experimental bounds. Fu¬ 
ture experimental searches for gamma lines are expected 
to have very good energy resolution, see for example the 
proposed GAMMA-400 [28], and one could hope to test 
this model in the foreseeable future. 


VI. SUMMARY 

In this letter we have shown that in simple models 
for Majorana dark matter one has a visible gamma line 
which could be used as a test of these dark matter models. 
In this context the dark matter is charged under a new 
gauge symmetry and the tree level annihilation into Stan¬ 
dard Model particles is velocity suppressed. The con¬ 
tribution to the gamma continuum spectrum from dark 
matter annihilation is suppressed and the gamma line is 
distinguishable from the continuum. 

Our numerical results are in agreement with the 
bounds on the relic density, the direct detection lim¬ 
its and bounds from missing energy and leptophobic Z' 
searches at the LHC. In models with gauged baryon num¬ 
ber and a self-conjugate dark matter candidate a striking 
line signature is a generic feature. Our novel results point 
toward an ideal scenario where using complementary ex¬ 
periments one can test a dark matter model. 
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